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Abstract: The spectral function of rho meson in hot /dense hadronic matter is 
studied by taking into account the nucleon-loop on quantum hadrondynamics model 
level. Different from the hot pion gas effect which changes the spectral function 
slightly, the nucleon-antinucleon polarization (Dirac sea) makes the spectral function 
very sharp and shifted towards the low invariant mass region significantly due to the 
decreasing effective nucleon mass. 



KEYWORDS: [Spectral function, vector meson, CjHD, finite temperature field theory 



It is widely prospected that in ultra-relativistic heavy-ion collisions a new phase 
of matter — quark-gluon plasma (QGP) may be generated and the spontaneously 
broken chiral symmetry may be restored. Among the proposed various signals which 
can detect the phase transitions in strongly interacting matter, the electromagnetic 
signals-dileptons and photons are considered to be the clearest ones, since they can 
penetrate the medium almost undisturbed. Considering the fact that light vector 
mesons can directly decay to dilepton pairs, the study about p meson in medium 
is especially interesting because the lifetime of p is smaller than that of the fireball 
formed in heavy- ion collisions and the decay of p can be completed inside it. The 
strong enhancement of low mass dileptons in central A — A collisions observed by 
CERES-NA45 has excited many theoretical works]!], @. Among them, the mass 
dropping of p mesons provides a good description of the data[||, |, || §]. 

The mass spectrum of hadrons in extreme hadronic environment has been widely 
discussed. One of the familiar results is the Brown-Rho scaling law indicating that 
the masses of nucleons and mesons decrease in hot /dense hadronic environment. 
For light vector mesons, it is difficult to extract an unique conclusion from different 
model calculations. The mean field theory (MFT) and effective Lagrangians can give 
similar results as the Brown-Rho scaling [|], [7]] , but the results from some QCD sum 
rules may be different from it and even contrast to it|| ^ |10[ . The property of p in 
hot/dense environment is not yet clear and needs further study. 

The hot topic in recent years is discussing the low invariant mass dilepton pro- 
duction in heavy ion collisions. The dilepton distribution is related with the spectral 
function of p meson, i.e., the imaginary part of the retarded propagator of p in 
medium [fLl , 12, l3| , which reflects the comprehensive medium effects. From the 



point of view of the partial chiral symmetry restoration the p spectral function has 
been widely discussed in recent literature such as in Refs.fi], |13| (and references 
therein). 

The key point of discussing the p spectral function is how to calculate the p 
meson self-energy in nuclear environment. In this paper, the effects of the decreasing 
effective nucleon mass on the p spectral function in medium are discussed. With 
quantum hadrodynamics model(QHD) of nuclear matter, the effective nucleon mass 
drops down with density/temperature and should be taken into account in discussing 
the effective mass as well as the decay width of p in medium, i.e., the spectral 
function. H. Shiomi and T. Hatsuda have found that the effective p meson mass 
decreases in nuclear matter due to the nucleon polarization at finite density with 
QHD-I type effective Lagrangian||14||. By extending their work to finite temperature, 



the temperature and density dependence of the p-meson spectral function due to 
nucleon loop is discussed. 

In Minkowski space, the self-energy of the light vector meson p can be expressed 
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where k 2 = k^ — k . The P£"(k) and P^ v {k) are the standard longitudinal and 
transverse projection tensors 

P°° = po, = p*0 = 0> p» = 5 ij _ W 

Jv 

pf = ^r- - aT - Pf- (2) 

The longitudinal and transverse elements of the self-energy can be determined by 

u L {k) = ^n 00 (fc), n T (*) = i/$ny(*). (3) 

Noting that the self-energy is related to the full (T>) and bare (T>q) propagators by 

W v = {V~ X Y V - (Vq 1 )^, (4) 

one can obtain 

pr pr 

j)^ v = k L (5) 

k 2 — m 2 — Hl k 2 — m 2 — Ut m 2 k 2 ' 

where m p is the free mass of p in vacuum. 

The property of in-medium p is determined by its full propagator. The imaginary 
part of the retarded propagator is referred to as the spectral function and related to 
the dilepton production with vector meson dominance model. The study of the p 
spectral function is attributed to calculating the in-medium self-energy. By concen- 
trating on decreasing nucleon mass effects, the in-medium p property is analysed in 
terms of the Feynman diagrams displayed in Fig.|l|. The pion effect has been ana- 
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Figure 1: The p self-energy Feynman diagrams associated with pnir and ppirir as well as 
pNN interactions. 



lyzed in Ref . [T1J , and in Minkowski space the pion- loop contribution to p self-energy 
can be written as 

UP n, (k) _ _ 2 T V- f d 3 p {2p + k),{2p + k) v , [Sp 1 

iv w v z^j ^Y(p 2 -mi)[(p + ky-miy g ^ gp ^ J i 



(27r) 3 p 2 — m\ ' 
(6) 



-2- 



where the O-component of n loop momentum is related to temperature via po = 
2mrTi, with n = 0, ±1, ±2, • • •. 

The contribution of nucleon excitations through nucleon-loop to p self-energy is 
analyzed in terms of the effective Lagrangian[14, 16] 

C p nn = g pNN (^T^V» - JZL-Va^rVf), (7) 

where V£ is the p meson field and ^ the nucleon field. It can be written as 



Po 



A' 



(8) 



where T M is defined as T^k) = 7^ + ^a^k u with cv = Mjy and M* N 

are the nucleon masses in vacuum and medium, respectively. The 0-component of 
nucleon loop momentum is related to temperature via po = (2n + l)iiTi + p* with 
p* being the effective chemical potential of nucleon. 

With residue theorem the pNN contribution can be separated into two parts 



ngTO = Fgr - 9M NN (k) + <>)• 0) 



The first part Il^ VAr (A;) corresponding to U F of Ref. [14| contains divergent integrals 
and needs to be regularized. However, it can not be renormalized by the conventional 
method. Using the dimensional regularization method and taking a phenomenologi- 



cal procedure [14, 17, 18, IS 



cPIlf N {k)/d{k 2 Y\ M ^ MN ^ =ml = 0, ( n = 0, 1, 2, • • • , 00 ) (10) 
one can eliminate the real vacuum (T = 0) part and obtain the finite result 

nf»(k) = e 9 -^- (p 1 + + ( _^_ )2 fc 2 A + A^ 2 p 2 \ ^ 



f,z,(t)W - n 2 \ AL ^ 2M N Z ' K 2M N _ / 

Pi— dxx(l — x) lnc, P2 — dxlnc, 
Jo Jo 

_ M* N 2 - x(l - x)fc 2 
C ~ M^ 2 -x)ife 2 " 

It is easy to see that the medium effect on U P J? N is also introduced by the effective 
nucleon mass in hadronic environment which will be determined by QHD-I 
below. 

The second part Il^f in @ is explicitly related to distribution function 
n^(T, p*) (corresponding to the Fermi sea contribution at T = 0). Its calculation is 
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lengthy. We list here only the longitudinal results 
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k 2 (A; 2 - 4p 2 ) + (k 2 - 2ook ) 2 , k 2 (P - 4p 2 ) + (fc 2 + 2a;fco) 2 



4p|k| 

where a, 6, u, n N , n N are defined by 
k 2 - 2p|k| - 2fc w 



In a + 



4p|k| 
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Before making further discussions, it should be emphasized that the inclusion of 
the pion tadpole diagram is very crucial for the current conservation condition(gauge 
invariance) [JTTJ] . It is interesting that the self-energy is divided into two parts such as 
in ([|) with residue theorem and the two parts separately satisfy the current conserva- 
tion. The current conservation condition ensures the transversality of IP", which is 
reflected in the properties of the projection tensors P£ v and P^ u such as k^P^^ = 0. 

After analytical continuation according to ko — > E + ie, E = \JM 2 + k 2 with 
M p being the p invariant mass in medium, one can obtain the real and imaginary 
parts of the self-energy. Then by using the definition of p meson spectral function 
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where Hl(t) is the total longitudinal or transverse self-energy, one can analyze the 
medium effects on p meson. 

The effective mass Mjy and effective chemical potential p* of the nucleon in ]JP NN 

~1. The QHD-I 



are determined from the simplest QHD-I at mean field level [20. 21 



is a renormalizable model including nucleon and meson degrees of freedom 



C=$b l *{id*-g v V> l )-(M N -g a 



^ + -(^0-m^) 



\f, v F^ + l -m 2 W + AC, 



(15) 
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where ip is the nucleon field, the neutral scalar field with free mass m s , = 
duVv—dvVfj, the field tensor of uj meson with free mass m v , and 5C the renormalization 
counterterm. 

In mean field approximation, the meson field operators are replaced by the 
ground-state expectation values 

cf) ^< (f) >= O) V» ->< V» >= g^V°. (16) 

With MFT by neglecting the vacuum fluctuation contribution, the classical field 4>q 
is a dynamical quantity which makes the nucleon mass drop and can be determined 



self-consistently bypl 22 



M N -M* N = g s <p = g 2 s < ^ >= ^-p s , (17) 



where the scalar baryon density p s is defined as 

7 f «,M* N 



Ps = J d'k-^- [n N (p*, T) + n N (ji*, T)] , (18) 

with the nucleon degenerate factor 7 = 4. 

The effective chemical potential fi* is introduced by 

P* = P~ 9vPN/ml, (19) 

with the baryon density 

P» = -fin* I d 3 ^[n N (p\T)-n N (p\T)}. (20) 
[2-kY J uj 

One can solve the coupled equations (|T7|) and (|19|) numerically. With the pa- 
rameters g* = 109.626, g 2 v = 190.431, m s = 520 MeV, M N = 938 MeV, and 
m v = 783 MeV according to the bulk binding energy and normal density of nuclear 
matter at saturation, the effective nucleon mass and chemical potential with MFT 
are displayed in Fig.^. 

The longitudinal spectral functions (|H ) in medium are drawn in Figs.|3] and f| for 
different hadronic environment. The parameters are chosen to be m n = 139.6 MeV, 
m p = 770 MeV, g 2 p = 4vr * 2.91, g 2 pNN = 6.91, k p = 6.1[|T§. It should be noted that, 
different from uj, the tensor coupling is very important for the interaction of p with 
nucleons. Here we use the coupling constants g P NN and k p determined from the fitting 
of the nucleon-nucleon interaction data done by the Bonn group. As shown in Fig.[3], 
when nucleon- loop is not taken into account, there is only temperature effect and 
the p meson spectral function changes weakly. With the increase of temperature, the 
spectral function becomes a little bit wider and is shifted towards the high invariant 



mass region slightly. This behavior agrees with the result of Ref. [11 . 
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Figure 2: Effective nucleon mass and chemical potential p,* determined by QHD-I 
as functions of scaled density (left panel ) and temperature (right panel ) with MFT. The 
normal nuclear density is chosen to be p$ = 0.16 /m -3 . 



The spectral function including the 
nucleon-loop contribution is shown in Fig 
The solid lines are still the calculation in 
vacuum, but the dashed and dot-dashed 
lines contain the density and finite tem- 
perature effects. The dotted line rep- 
resents the result of pure density effect 
(T = 0, pb 7^ 0). The density effect on p 
spectral function is obviously important 
by comparing Fig.|3| with Fig.^. The nu- 
cleon excitation shifts the spectral func- 
tion towards low invariant mass region. 
As indicated in the left panel of Fig||, in 
the case of low density, Pn/po — 0.1, the 
density effect is already remarkable com- 
pared with the pure temperature effect 
shown in Fig.|3|. The peak position of spectral function is shifted from 0.77 GeV in 
vacuum to about 0.63 GeV in medium. 

The second characteristic of nucleon loop effect is the sharping of p spectral 
function. Even for low density situation, see the left panel of Fig|| with pn/po = 0.1, 
the p meson width is already reduced to about half of that in vacuum. When density 
is high enough, see the right panel of Fig.f| with Pn/po = 1-5, the spectral function 
approaches to a 5-like function. This is mainly a consequence of the reduced two-pion 
phase space due to the p mass decrease in medium. 

Another significance of the nucleon loop is its approximate saturation at high 
density. At Pn/po = 1-5, the two lines with temperature T = 0.1 GeV and 0.15 GeV 
almost coincide with each other. When density goes up further, for instance taking 
Pn/po — 2, the spectral function is almost the same as the corresponding line in the 
right panel of Fig.^ for Pn/po — 1-5. This saturation tendency can be seen from the 
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Figure 3: The p meson spectral function 
without considering the nucleon-loop at mo- 
mentum k = 0.75 GeV. The solid line rep- 
resents the vacuum situation, while dashed 
and dot-dashed lines correspond to temper- 
ature T = 0.1 and 0.15 GeV, respectively. 
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Figure 4: Spectral function including the nucleon-loop contribution at momentum k = 
0.75 GeV. Line style is similar to Fig.|| except that the dotted lines represent the T = 
situation. Left panel: pn/po = 0.1; Right panel: pn/po = 1-5. 
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Figure 5: Effective p meson mass vs scaled baryon density for T = 0.15 GeV with MFT 
(left panel ) and RHA (right panel). For comparison, the effective nucleon mass line is also 
drawn out. The dashed lines correspond to T = situation, while the solid lines to those 
in the hot medium. The symbol "1" represents the case considering form factor, while "2" 
without the form factor contribution. 



effective mass m* determined by the pole position of its propagator in medium |[L4 



k 2 Q - m 2 - lim U L{T \k) = 0, (21) 



p Ikl— o 



as indicated by Fig|| For comparison, we also give the effective nucleon mass Mjy. 
By taking MFT, m* and are displayed in the left panel of Fig. |^. When density 
is higher than po, m* approaches its minimum value and then increases with the 
increase of density a little bit. The result of m* with the nucleon propagator obtained 
with the relativistic Hartree approximation (RHA) by including vacuum fluctuation 
contribution instead of MFT is quite similar except the numerical difference, as 
indicated in the right panel of Fig|| Vector meson masses will decrease with the 
increase of density for a model including the polarization of nucieon-antinucleon with 



a decreasing effective nucleon mass[14, 21|. The numerical calculation indicates that 



the Dirac sea contribution dominates. 

It should be noted that the interactions between mesons and nucleons are point- 
like in QHD. Considering the fact that mesons and nucleon are composite, a phe- 
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nomenological monopole form factor is multiplied at each meson-N vertex 



F a {k>) = (22) 

where is the four-momentum transfer, m a the mass of the exchanged meson and 
A a the relevant cut-off mass. For pNN vertex, A p is usually taken as A p = 1.5 GeV. 

The inclusion of form factor attenuates the behavior of effective mass m* sig- 
nificantly, which makes the magnitude of mass shift in medium smaller than that 
without the form factor. For analyzing the influence of form factor on the effective 
mass m*, the two situations with/without form factor contribution have been drawn 
out in Fig.||. The displayed spectral function in Fig. |] has included the form factor 
contribution. 

Our result of the spectral function is different from previous results such as in 
Refs.[|], [T3|]. From Figs.|] and [5], one can see that effects of the decreasing effective 
nucleon mass on p spectral function are important but they are not considered in 
Refs.[||, [13| . On the other hand, by focusing on the effects of decreasing nucleon mass 
on QHD-I level, any other decay channels which can contribute to the imaginary as 
well as the real parts of the self-energy are not included in our calculation and also 
cause the difference. 

In summary, we discussed the density and temperature effects on p meson spec- 
tral function with the effective Lagrangian. Different from the pure temperature ef- 
fect of the pion loop which changes p meson mass only a little and shifts the spectral 
function towards high invariant mass region slightly, the nucleon-loop changes the 
p meson properties differently and significantly: 1) The spectral function is shifted 
towards low mass region remarkably even in the case of low density; 2) Together with 
the shift towards low mass region, the spectral function becomes very sharp due to 
Dirac sea polarization with the decreasing nucleon mass in medium; The density 
effect on the spectral function approaches to saturation when the density increases. 
3) The influence of a phenomenological vertex form factor may prevent the effective 
p mass decreasing to some extent within hadronic level. Although the behavior of 
m* in low density region is expected to be not so stiff, our result indicates that the 
effects of decreasing nucleon mass may be remarkable and should be taken into ac- 
count in discussing the in-medium p property. The relation between our result and 
those from QCD sum rules esp. in the low density region is not yet well understood 
and needs to be clarified further. 
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